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Non-equilibrium quasi-condensates in reduced dimensions 
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Abstract -We develop a generic theory of the fluctuations on top of a non-equilibrium Bose- 
Einstein condensate in reduced dimensions. Analytic expressions for the momentum distribution 
of the non-condensed cloud and for the long-distance behavior of the spatial coherence are ob- 
tained for the different cases. Comparison of our predictions with available experimental data on 
condensates of exciton-polaritons and on surface-emitting planar laser devices is finally made. 
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Non-equilibrium phase transitions are among the most 
exciting topics of non-equilibrium statistical mechanics [I] . 
In the last decades they have been studied in a variety of 
different physical systems, of either classical ,2 -4 or quan- 
tum 0[6] nature. However, as compared to their equilib- 
rium counterparts, much less is known about their gen- 
eral features, in particular for what concerns the critical 
behavior in the vicinity of the transition point. A most 
important class of non-equilibrium phase transitions take 
place in systems at the so-called non-equilibrium steady 
state (NESS) where a dynamical balance of driving and 
dissipation replaces the usual thermal equilibrium condi- 
tion of standard equilibrium statistical mechanics. 

Even though laser operation has played a central role in 
most developments of contemporary experimental atomic, 
molecular and optical physics, its potential as a workbench 
for non-equilibrium statistical mechanics studies has not 
been fully exploited yet. The interpretation of the laser 
threshold as a second-order phase transition dates back 
to the early 1970's with the pioneering works by DeGior- 
gio and Scully [7] and by Graham and Haken [8j|9]: the 
order parameter of the transition is the amplitude E of 
the electromagnetic field in the laser cavity, which gets 
a well defined phase by spontaneously breaking the U(l) 
symmetry corresponding to phase rotations, E — >• E e tS . 

While most textbook discuss this analogy in the case 
of single-mode cavities only, the spontaneous symmetry 
breaking phenomenon is -rigorously speaking- restricted 
to spatially extended systems. Only in this case, it is in 



fact meaningful to take the long-distance limit involved in 
the Penrose-Onsager definition of an ordered state, 



lim (#(x)£(x')>^0. 

|x— x' |— >oo 



(1) 



Most remarkable examples of spatially extended laser de- 
vices are the so-called vertical cavity surface emitting 
lasers (VCSELs), whose planar geometry makes a con- 
tinuum of in-plane modes to be available to both lasing 
and fluctuations. Even though such VCSEL devices have 
attracted a strong interest from the point of view of non- 
linear optics and of all-optical information processing |10| , 
very little work has been devoted to the non-equilibrium 
statistical mechanics features of their operation. On the 
other hand, a great attention has been devoted in the 
last years to studies of the closely related Bose-Einstein 
condensation (BEC) phase transition in gases of exciton- 
polaritons in semiconductor microcavities In the 
strong light-matter coupling regime, the cavity photon 
is strongly mixed to an excitonic transition in the cavity 
medium, giving rise to new long-lived bosonic excitation 
modes, the so-called exciton-polaritons: for strong enough 
pump intensities, the polariton density reaches the thresh- 
old for BEC and long-range spatial and temporal coher- 
ence appears. Recently, BEC was also observed in a pho- 
ton gas in a macroscopic cavity filled of dye molecules [12] . 

Even though most features of these BEC experiments 
are accurately described under the assumption that ther- 
mal equilibrium is reached within the photon/polariton 
gas on a time-scale shorter than the particle decay rate 



p-1 



A. Chiocchetta et al. 



7, several authors have pointed out intriguing conse- 
quences of the non-equilibrium nature of these optical 
BECs, whose steady-state is determined by a dynami- 
cal balance of pumping and particle losses: the Goldstonc 
mode associated to the spontaneously broken U(l) sym- 
metry has a diffusive rather than sonic nature [T5HT5] and 
the momentum-space shape of the condensate strongly de- 
pends on the pump spot geometry [T6iri7] . 

As the spontaneously broken symmetry is the same U(l) 
in all cases of polariton BEC, photon BEC, and laser op- 
eration, it is possible to construct a unified theoretical de- 
scription of these phenomena under the general name of 
non- equilibrium Bose- Einstein condensation: the role of 
stimulated emission of photons by a population-inverted 
laser medium is played in the BEC case by stimulated 
scattering from the thermal cloud into the condensate 
mode; the actual distance from an equilibrium state may 
of course depend on the microscopic details of the system. 
The purpose of this Letter is to develop such a general 
theory and then to apply it to the long-range coherence 
properties of (quasi)-condensates in reduced dimensions: 
independently of collisional thermalization processes, we 
shall see that fluctuations in the non-equilibrium steady- 
state recover the ones of the corresponding equilibrium 
system at a temperature determined by the characteristic 
decay rate. 

Stochastic Gross-Pitaevskii Equation. Among 
the different theoretical frameworks have been used to 
study non-equilibrium BEC transitions in planar optical 
devices [HJ[T51[TII] ; m this Letter we shall use the simple, 
yet powerful techniques based on the so-called stochastic 
Gross-Pitaevskii equation (SGPE). The idea is to start 
from the mean-field equation for the coherent condensate 
field and then to add a stochastic noise term to include 
quantum and/or classical fluctuations. For the case of 
a coherent pump, this procedure can be made rigorous 
within the Wigner representation of the quantum field and 
has been successfully applied in Monte- Carlo studies of 
the field coherence across the transition point [T3l[20] . A 
first successful extension of this method to the incoher- 
ent pump case that is under investigation here was made 
in [21j on phenomenological grounds: a more rigorous 
derivation requires a quantum description of the ampli- 
fying mechanism that will be published in a future, more 
extensive work based on [22] . 

Independently of the microscopic nature of the Bose 
field </)(x, t) (which can describe equally well the electro- 
magnetic field of a VCSEL and the polariton field of a 
polariton condensate) , its low-energy physics can be accu- 
rately described by a d-dimensional SGPE of the form [T5] : 
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where m is the boson mass of the field, 7 is the damp- 
ing rate, Pq is the pumping rate and n s is the saturation 



density of the pumping mechanism. As usual, the validity 
of the SGPE model requires that the nonlinear coupling- 
strength g > describing repulsive contact interactions 
between the particles be weak enough for the gas to be 
dilute. 

As usual in open systems, the presence of dissipative 
terms in the motion equations requires including corre- 
sponding stochastic terms dW(x.,i). For the sake of sim- 
plicity, we phenomcnologically choose the simplest Ito 
form of a temporally white complex Gaussian noise term 
with a random phase and spatially local correlations, 

{dW(x,t)dW*(x',t)) = 2A w ,5 (d) (x-x')*, (3) 
(dW(x,t)dW(x.',t)) = 0. (4) 

The assumption of a white, temporally uncorrclated noise 
is based on the physical assumption that the characteristic 
time of the bath is much shorter than the characteristic 
time of the system. A brief account of effects beyond the 
white bath approximation will be given in the last section 
of the Letter. An ab initio calculation of the diffusion 
coefficient will be presented in a future publication: 
based on previous work on the OPO pump case |20j . we 
expect a value « 7 for pump strengths ranging from 
zero to not too much above the transition point. 

Mean-field approximation and collective excita- 
tions. — The first step to understand the behavior of 
the system is to perform a mean-field approximation and 
look for the stationary state of the deterministic GPE 
that is obtained by neglecting the noise terms in ([2]). For 
simplicity, we restrict our attention to the case of a spa- 
tially homogeneous system, for which we can look for a 
translationally symmetric mean field solution of the form 
0(x,t) =<t> e- iut . 

Independently of the value of the pump intensity Po, 
the trivial field <j> = is always a solution of the GPE. For 
low pump intensities below the threshold Po < 7: this is 
the only solution and is dynamically stable. On the other 
hand, for pump intensities P > 7 above the threshold, it 
becomes dynamically unstable and is replaced by another 
solution with a finite density |0o| 2 = ^s(P) — 7)/7 & n d an 
oscillating frequency u = u>o + g\4>o\ 2 '- the global phase of 
4>o acquires a well-defined value by spontaneously breaking 
the U(l) symmetry of J2j) . In the analogy with second 
order phase transitions, the order parameter hallmarking 
the transition is then the complex field amplitude </>o- 

The linearization of (J2J around the finite solution de- 
termines a pair of coupled Bogoliubov partial differen- 
tial equations for the fluctuation field <5</>(x, t) defined as 
(/>(*, t) = [0o + t)] e -1 "*, and its complex conjugate 
6<fr*(x,t). These fluctuations describe the excitations of 
the Bose field on top of a pure condensate solution, that 
is the non-condensed fraction of the gas. 

In the spatially homogeneous case under investigation 
here, it is useful to work in Fourier space and separate the 
evolution of the different in-plane components. For each 
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Fig. 1: Bogoliubov dispersion of collective excitations of the 
non-equilibrium condensate as a function of the (real) in-plane 
momentum k. Upper panels show the real part of the fre- 
quency, lower panels show the imaginary part. The pan- 
els (a-b) in the left column are for the case of a frequency- 
independent bath. The panels (c-d) in the right column model 
the frequency-dependence of the bath as discussed in the last 
part of the Letter. System parameters: Po = 27, gn a = 7, 
Uc = 27, LOQ = 0. 



(real-valued) k we are then left with a pair of coupled 
linear differential equation for <$0(k) and [5<j){— k)]*, 



W(-k))* 



ty(k) 

dt v w-k)r j — v w(-k)r 

with the k-dependent Bogoliubov matrix Ck defined as 
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in terms of the free particle dispersion ek = hk 2 /2m and 
the interaction energy /1 = g \ <fio\ 2 - By simple diagonaliza- 
tion of £k, one gets a double-branched excitation spectrum 
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where = ^/ £ k(ek + 2/i) is the standard Bogoliubov dis- 
persion of equilibrium systems and the effective damping 
rate T = ^(Pq— 7) /Po. An example of this non-equilibrium 
Bogoliubov dispersion is shown in the left panels of FigQ] 
as first discussed in [T3HT5] its most remarkable feature 
is the diffusive nature of the Goldstone mode describing 
long- wavelength phase twists in the condensate: the real 
part of its frequency vanishes within a finite region of 
in-plane wave vectors around k = 0, while its imaginary 
part grows quadratically from 0, u;^ ~ —ih/j,k 2 /mT. The 
other branch ui^ describing density fluctuations (intensity 
fluctuations in the language of laser physics) is instead 
gapped and tends to — iT for low k: remarkably, the ef- 
fective damping rate T vanishes when the critical point is 
approached from above, while it tends to the bare damp- 
ing rate 7 far above the threshold. 



Bogoliubov theory of linearized fluctuations. 

The linearized field equations ^ are the starting point to 
study the fluctuations of the Bose field around the pure 
condensate solution. As usual, this approach is accurate 
when fluctuations are small, that is far away from the 
critical point. A more complete description of the large 
fluctuations in the critical region requires more sophisti- 
cated approaches and will be postponed to further work. 
A first remarkable work in this direction making use of 
rcnormalization group techniques appeared very recently 
in [23]. 

Including back the noise term into their RHS, one ob- 
tains for each k a pair of linear Ito stochastic differential 
equations from which it is straightforward to compute all 
the steady-state correlations of the field, e.g. the momen- 
tum distribution of the particles 



nft = (2 7 r) d |0 o | 2 ( 5W(k) 
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in experiments, this quantity is directly accessible as the 
angular distribution of the far field emission. The former 
term in (J5J describes the condensate, while the latter one 
account for the non-condensed cloud. Note that the UV 
divergence due to the second term in the bracket (the one 
with no k dependence) is an artifact of the white bath 
assumption and disappears as soon as the frequency de- 
pendence of the bath in included in the theory 

In what follows, we shall instead concentrate on the low- 
k, long-distance physics, where our effective theory is ac- 
curate. The low-k behavior of n£ s in the different cases 
is summarized in the table and compared to the predic- 
tion of the standard Bogoliubov theory [21] for equilibrium 
systems at T = and at T > 0. 
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In the presence of interactions g > 0, the behavior of 
the non-equilibrium system recovers the one of the corre- 
sponding equilibrium system at an effective temperature 



T, 



eff 



D^j^+T 2 ) 
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(9) 



independently of the absence of thermalizing collisions be- 
tween different Bogoliubov modes, the momentum distri- 
bution recovers the typical shape of a gas at thermal equi- 
librium, n k s oc k~ 2 . As a result, the mere observation of 
a thermal-like momentum distribution does not appear to 
be an unambiguous proof of a thermal equilibrium state 
in the polariton gas. The situation is very different in 
the case of a non-interacting gas with g = 0, where the 
momentum distribution scales as n£ s oc k~ A . 



From the momentum distribution 



l k 1 



it is immediate 



to obtain by Fourier transform the one-time spatial corre- 
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lation function 



G(|x- 



<0*W^(y)> 



d d k 

(2n) d 



a ik(x-y) 



(10) 



Long-range order is then denned according to the Pcnrose- 
Onsager criterion by the condition 



lim 
|x-yK 



G(|x-y|) 
G(0) 



^0. 



(11) 



While the UV divergence in (fTU|) is an artifact of our 
white bath model and can be safely removed by means 
of a UV cut-off, the convergence vs. divergence of the 
integral on the IR side determines the presence vs. ab- 
sence of long-range order. Inserting the explicit formula 
(J8j> into (fT0|) . one immediately sees that the long-range 
order in an interacting gas is only stable in d > 3, while 
an even higher d > 5 is required in a non-interacting gas. 
Once again, the prediction for the interacting gas recovers 
known results for the equilibrium gas and can be seen as 
a non-equilibrium version of the Mcrmin- Wagner theorem 
of equilibrium statistical mechanics [25| . 

Density-phase Bogoliubov theory. In \ow-d 
where the condensate is destabilized by fluctuations, the 
Bogoliubov theory sketched in the previous section is in- 
consistent in the thermodynamic limit, as it implicitly as- 
sumes the existence of a condensate to start with. To go 
beyond this difficulty, more sophisticated versions of the 
Bogoliubov theory have been developed, based on the idea 
of a quasi- condensate [261127] : The basic idea is to separate 
fluctuations in density and phase fluctuations according to 



(x,t) = ^n + Sn(x,t) e ie(x '* } e~ iut . 



(12) 



We then assume that density fluctuations Sn around the 
average value hq are weak Sn <C hq and that the phase 
9 smoothly varies in space. As the phase 9 can take ar- 
bitrarily large values, no assumption on the presence of 
a condensate is made and the theory is able to describe 
configurations with no long-range order. 

Within this density-phase Bogoliubov approach, the dif- 
ferent k components of density and phase fluctuations de- 
couple and evolve according to linear differential equations 
very similar to (j6|) including Ito noise terms. Averaging 
over noise, we obtain the following forms for the correla- 
tions of density and phase fluctuations at a given k, 



(Sn k Sn^) ss 
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where S kk > = (2Tr) d S < - d) (k - k'). 



(13) 

(14) 
(15) 



Spatial coherence in reduced d. — Fourier trans- 
form of (fTU|) back to real space shows that the correlation 
function of density fluctuations quickly tends to zero on a 
length scale set by the usual healing length £ = \JK 2 / /im. 
At larger distances, the field correlation functions can then 
be written in terms of phase fluctuations only. As these 
inherit the Gaussian statistics of the noise, we can rewrite 
the field correlation as 



G(x - y) ~ n e 



-G„ e (|x-y|) 



in terms of the phase correlation function 
G ee (x-y) = i([0(x)-0(y)] : 



(16) 



(17) 



Focussing on the long-distance behavior of the spatial 
correlation function, we can safely neglect the short range 
correlations due the first term in Q14[l and write the cor- 
relation function in the form 



d d k e lk x - 1 



(2w) d 



k 2 

d d k e lk x - 1 



(27r) d k 2 + £- 2 
where the additional length scale £o is defined as 



n T 



(18) 



(19) 



In the following of this Section, we shall separately con- 
sider the physically relevant d = 1,2,3 cases, and highlight 
the peculiarities of each of them. 

d = 3. As we have already mentioned, in d = 3 both 
Gee and G correlation functions have finite long-distance 
limits as soon as the interaction strength is finite g > 0: 
this signals the existence of long-range order, while the 
momentum distribution of the non-condensed particles fol- 
lows the typical k~ 2 law of an equilibrium system at the ef- 
fective temperature T e fj determined by ([9]). On the other 
hand, in the opposite case of a non-interacting gas with 
g = 0, taking the \i — > limit leads to a field correlation 
function 

G{x)(xe~^^ (20) 

that exponentially decays at large distances within the 
characteristic length £o- This is, of course, in stark con- 
trast with the existence of a BEC transition in the text- 
book case of a d = 3 ideal Bose gas. 

d = 2. In this case, the long-distance behavior of the 
correlation function can be evaluated making use of known 
properties of the Bessel functions to have the power-law 
form 



G(|x 
with an exponent 



const - 



1 £ 2 
2^Co 2 (M) 
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(21) 



(22) 
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This power law decay behaviour is characteristic of 
a superfluid (but not Bose-Einstein condensed) two- 
dimensional Bose gas at temperatures below the 
Bcrczinskii-Kosterlitz-Thouless transition [55] and was 
first discussed in a non-equilibrium context in [14] . Iden- 
tifying the superfluid density n s with the average den- 
sity no and the temperature T with the effective temper- 
ature T e ff defined in (J9)), the form ([22]) of the exponent 
77 matches the usual form tjekt = 1/psAj., for equilib- 
rium BKT gases, with the de Broglie wavelength defined 
as usual as At = [2irh 2 / inksT] 1 / 2 . In spite of this, it is 
essential to note that, while in the equilibrium case ther- 
modynamic stability of the quasi-ordered BKT phase im- 
poses an upper bound 77 < 1/4 to the exponent [28], no 
such condition is known for the non-equilibrium case and 
larger exponents might be observable. Experimental in- 
dications in this direction were recently reported for an 
interacting polariton gas in [29| . 

On the other hand, the field correlation function of a 
non-interacting g — > gas has the form 



Gn 



l-£ + ln2-ln( - 



£o(0) 



(23) 



which quadratically grows in the r 3> £0 limit. Q is here 
the Euler-Mascheroni constant, Q ~ 0.577; the logarith- 
mic correction remains finite as the g = condition cor- 
responds to a diverging value of the healing length £. As 
a result, all the quasi-order that was present in the BKT 
state disappears in the g —> limit and is replaced by 
a much faster almost Gaussian decay: experimental re- 
sults shown in the same article [29] for the very weakly 
interacting photon regime suggest that in this case the 
correlation function indeed decays in the quasi- Gaussian 
way predicted by ([23]) . 

d = 1. Straightforward integration of (JTSJ) in the d = 1 
case predicts a fast, diverging growth of Gee of the form 
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.(*) = \ 



MJ t Uo(o) 
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1) 



(24) 

Depending on the relative value of the two characteristic 
distances £ and £qj this leads to different functional forms 
of the long-distance tail of G(x). 

For £0 ^ £ (which is the case of an interacting gas at 
strong pump values well above threshold), the field corre- 
lation function has an exponential tail G{x) ~ q-\ x \/^ as 
first predicted in [13] for the OPO pumping scheme. This 
exponential form closely follow the corresponding case of 
a d = 1 gas at equilibrium at a finite T > [30] . Also 
the characteristic length t\ has the same form in the two 
cases: the explicit expression 



2no/ir 



m%(/i 2 +r 2 ) 



(25) 



matches the equilibrium one IjP = nX^/ir upon identify- 
ing the temperature T with the effective temperature T e g 



defined in The situation is somehow different in a 

non-interacting gas and, more in general, for £0 £: in 
this case, the long-distance tail can be approximated with 
a Gaussian form G(x) ~ e - ^^ 2 ' of characteristic length 
^2 = 2\/£o(0)/£- Experimental indications of a reduced 
coherence in one-dimensional polariton condensates was 
reported in [3"Tll32j . 

Frequency dependent pumping. All the discus- 
sion so far has been based on a model assuming amplifica- 
tion and dissipation baths with no frequency dependence. 
This simplifying assumption was based on the physical in- 
tuition that the correlation time of the bath should be 
very short, which justifies a Markovian approach based on 
Ito partial differental equations. More sophisticated ki- 
netic descriptions of the amplification mechanism [33,34 
as well as experimental observations of condensation in 
harmonic traps [35] suggest the presence of some mecha- 
nism tending to concentrate the emission into low-energy 
modes and quench the high-energy ones. 

Following [36], a simple way to model this behavior is 
to include some frequency-dependence of the amplifica- 
tion on a characteristic scale u c , e.g. by replacing the 
amplification coefficient Pq with the differential operator 
P = Pq [1 + cj f T 1 (u;o — id/dt)]: in this way, amplification 
(and then condensation) into low energy modes is favored, 
while it is completely suppressed above the cut-off fre- 
quency uj > uj c + ujq. As a result, the transition point 
is still at Pq = 7, but the competition of the frequency- 
dependent amplification with the interaction-induced blue 
shift makes the condensate density to saturate at large 
pump powers, \(j)\ 2 = n s (P - PS)/{Pq + Pagn s /bj L ). 

An example of Bogoliubov excitation spectrum around 
the spatially homogeneous, plane- wave condensate solu- 
tion is shown in the right panels of FigQ] while the Gold- 
stone mode with a diffusive dispersion at low k's is still 
visible, the damping rate of the Bogoliubov modes rapidly 
grows at large k's. The momentum distribution of the 
excitations is obtained by adding phenomenological noise 
terms to the Bogoliubov equations as done in the first part 
of this Letter. This leads to the closed expression: 



r + M(e k + /i) 



(26) 



in terms of the parameters f = 7(Po — Pq )/[Po(l+<?n s /w c )] 
and M = P /[o; c (l + \(j) \ 2 /n s )}. While the low-k part 



of the momentum distribution (|26|) retains the same be- 
havior as in ((SJ), the frequency-dependent amplification is 
able to regularize the UV behavior of the theory, mak- 
ing all relevant integrals to nicely converge at large k. Of 
course a quantitative account of the details of the momen- 
tum distribution requires a more sophisticated model of 
the frequency-dependence, still our model is able to qual- 
itatively recover its fast decrease at large k as typically 
observed in Bose-Einstein condensation experiments. In 
particular, note how our theory predicts that this behav- 
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ior is completely independent of the collisional thermal- 
ization rate: this may explain the unexpected experimen- 
tal observation of a thermal-like distribution in a gas of 
non-interacting photons in a VCSEL device in the weak- 
coupling regime [37] . The role of this effect in determining 
the thermal photon distribution observed in the photon 
BEC experiment of [12] is under investigation. 

Conclusions. — In this Letter we have developed a 
fully analytical theory of the coherence properties of non- 
equilibrium Bosc-Einstein quasi-condensates in reduced 
dimensionalities. Even though the non-equilibrium sta- 
tionary state is determined by a dynamical balance of 
pumping and dissipation rather than by a standard ther- 
modynamical equilibrium condition, the momentum dis- 
tribution and the long-distance spatial coherence of a 
weakly interacting gas is qualitatively similar to the one of 
the corresponding equilibrium system at a finite temper- 
ature. Most remarkably, this conclusion does not depend 
on the collisional thermalization rate in the gas and, with 
some modifications, extends to the case of non-interacting 
particles. Thanks to the generic nature of our theory, its 
conclusions can be applied to a broad range of systems 
displaying a non-equilibrium BEC-like phase transition, 
from laser operation in VCSEL devices, to photon BECs, 
to polariton BECs. As a next step, we will derive ex- 
plicit predictions for the amplification and noise parame- 
ters starting from microscopic models of the most relevant 
physical systems, so to provide quantitative predictions for 
the emission spectra and the coherence functions in the 
different cases. 
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